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The scaling properties of the maximal height of a growing 
self-afBne surface with a lateral extent L are considered. In 
the late-time regime its value measured relative to the evolv- 
ing average height scales like the roughness: h*j^ ~ L°' . For 
large values its distribution obeys log P{h}^) ~ — 
charaterized by the exponential-tail exponent a. In the 
early-time regime where the roughness grows as t'^ , we find 
hi ~ t^[lnL- (f ) Int -I- C]^^'' where either & = a or 6 is the 
corresponding exponent of the velocity distribution. These 
properties are derived from scaling and extreme-values argu- 
ments. They are corroborated by numerical simulations and 
supported by exact results for surfaces in ID with the asymp- 
totic behavior of a Brownian path. 



Kinetic roughening of self-affine surfaces has contin- 
uously attracted much interest in the last two decades. 
It is observed in many systems such as crystal growth, 
vapor deposition, molecular-beam epitaxy, electrochem- 
ical processes, bacterial growth, burning fronts, etc. |^] 
[H I§1 Its investigations mostly dealt with the scal- 
ing properties of the surface roughness. Other charac- 
teristics, which have been studied more recently, include: 
The distribution of the width ||^ |6j, the distribution of 
the site-average height velocity in the ID Kardar-Parisi- 
Zhang (KPZ) model [||, the density of the extrcma 
[H, the persistence and cycHng effects In the 

present work we focus on another important property: 
The scaling behavior of the highest point of the surface. 
In contrast with all characteristics listed above, the max- 
imal height is not a site-average property. Rather, it is 
an extremal one which takes just one single value for ev- 
ery realization of the surface. Since the maximal height 
is one of the two most uncommon values of the height 
of the whole surface (the second being the lowest height 
which follows a similar behavior), to obtain its scaling 
behavior and that of its distribution, an inherently dif- 
ferent approach is required. It blends together the ususal 
scaling methodology to the roughness with extreme- value 
statistics of rare events [|2| . 

Besides the theoretical challenge involved in rare events 
statistics of variables with intricate spatial and temporal 
properties, the maximal height is of a paramount im- 
portance in technological applications. It is the most 
significant property for systems with corroded surfaces 
in which the deepest, or the weakest, point determines 
the onset of breakdown |13]. Another such example are 



batteries in which a shortcircuit occurs when the highest 
point of the metal surface accumulated on one electrode 
reaches the opposite one. 

In addition to the absolute height measured from the 
substrate, we define the relative height as the height 
measured from the evolving average height of the sur- 
face. The scaling properties of the maximal height (MH) 
and the maximal relative height (MRH) are derived from 
the combined scaling and extreme statistics approaches. 
They are supported by numerical simulations and by rig- 
orous results from extremal excursion properties of Brow- 
nian paths. Our main findings are: In the late-time 
regime the MRH scales like the roughness and different 
properties of its distribution are elucidated. In the short- 
time regime the maximal height grows logarithmically 
slower than the roughness and has a logarithmic depen- 
dence on L (while the roughness is L-independent). The 
exact form of these logarithmic terms depends on the 
late-time distribution of either the velocity or the MRH. 

The surface roughness w{L,t), where L is the lateral 
size of the system is defined as w'^{L,t) = H{r,tY — 

2 

HL{t) . In this definition H{r,t) is the surface height. 
An overbar denotes average over the sites r of the 
basal plane of size (D — d — 1). In particular, 
HL{t) = L^^^P(zLD H{f,t) is the average height of a 
given surface at time t, and HL(t) = vt (with v being 
the growth velocity). We also define the roughness as 
W{L,t) =< w^{L,t) >^/^, where angular brackets de- 
note average over the ensemble of surface configurations. 
We recall that that it obeys the following scaling form 0| 
|l|: W{L,t) ~ L"£f(L/^(t)), where ^(i) - i^/^ is the lat- 
eral correlation length (z is the dynamic exponent). For 
large time t > tx (L) — U- it reaches its time-independent 
value: Wl ~ i", while for t < {L) it is L-independent: 
W{t) ~ i'^, where [3 — ^ is the growth exponent. Some 
generic models of surface growth are: Random depo- 
sition (RD), Edwards- Wilkinson (EW), KPZ 0, 
DasSarma-Tamborenea (DT) iQ etc. 
The relative height is given by: 

hL{r,t)^HL{r,t)-lhM (1) 

The height at a site is the sum of the average height 
and the relative height. If either one of them fluctuates 
much more than the other, it will make the dominant 
contribution to the spread of the height distribution in 
general and to that of the MH (denoted by H^{t)) in 
particular. 
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A most important distinction with regard to the MH 
is, therefore, between a first category of systems in which 
the velocity v at the large t limit, has a non-trivial distri- 
bution P{v) (because the instanteous rate of deposition 
on the surface is affected by its configuration), and a sec- 
ond category of systems in which P{v) is trivial because 
the instantenous deposition rate is indepnedent of the 
surface configuration. Typical systems of the first cate- 
gory are those in the KPZ universality class Q. In the 
second category are, e.g., systems in the EW and DT uni- 
versality classes with either a constant or a white-noise 
deposition rate. 

The rational for this categorization will become ap- 
parent in the discussion of the early-time regime which 
is addressed below following the the analysis of late-time 
regime. In this latter regime, systems from both cate- 
gories with P{v) ^ 5{v — vo) have their MH distribution 
dominated by that of the average height itself. Thus their 
height distributions are given essentially by those of their 
respective velocities (multiplied by the time). For such 
systems the MRH h*{r,t) = H*{r,t) - IhJf) describes 
the subdominant fluctuations and is of interest as well. 
The MH distribution for some important 2"'' category 
systems in which v is not fluctuating {e.g., electrodeposi- 
tion at a constant current) is obviously given by that of 
the MRH, up to a uniform shift by vt. 

We begin by analyzing the behavior of the MRH in the 
late-time regime. In this regime the probability distri- 
bution of the relative height is time- independent. Con- 
sistent with the scaling description, the roughness 
is the only relevant scale of this distribution. Thus the 
height distribution depends only on m = and P{u) is 
independent of L (by P{x) throughout we denote the dis- 
tribution of X, rather than one particular distribution). 
Therefore the MRH distribution should only depend on 
u* ~ h*^/wL, and be independent of L. This is confirmed 
numerically for the distributions of the ID EW model for 
different L, which are depicetd in Fig. 1. 

Our simulations indicate that this distribution van- 
ishes for h*^ 0, faster than exponential. For large 
arguments h'^ oo, the distribution also drops faster 
than exponential and could be fitted to a behavior ~ 
exp[{-A{hl/L°')°'}] with the exponent a = 2. Therefore 
all moments < (/i^/W^l)*' > (with k positive or negative) 
are finite. In particular, both the average and the typical 
MRH scale like Wl. 

While the MRH distribution for this model depends 
on the boundary conditions (BC), its functional de- 
cay at large values does not. We find a behaior ^ 
exp[{-A{hl/L°')''}] with a = 2 independent of the BC 
p6[ . (in contrast, the amplitude A varies with the BC). 
This is not the case for the behavior of P{h*j^/WL) near 
the origin which is BC sensitive. 

These numerical observations are reinforced by ana- 
lytical results for the ID Brownian path which faithfully 
describes the ID EW and KPZ surfaces in the late-time 



regime. We represent them as a Brownian path with fixed 
BC -ff(O) = H{L) {i.e. as a Brownian bridge), allowing 
the average height to fluctuate. For the statistical prop- 
erties of h{r) — H{r) — Hj^ of this Markovian process, 
these BC are equivalent to periodic ones. 

The statistical behavior of the maximal height may 
be reflected either by the MRH h\ = Hl~ H or by 
Ai = {H^ — Hy. More complete results are obtained 
for the distribution of A^. However, h*^ and A^ are 
simply related by: /i^^ — Al — w\. 

Averaging this equation over different surface realiza- 
tions yields: 



< h*' >=< Al > -Wl 



(2) 



Using path-decomposition techniques ||T^ , we flnd the 
Laplace transform (generating function) of the distribu- 
tion P(Al): 



P{\) =< exp{-\AL} >= 



V2XL 



3/2 



(3) 



This immediately implies that LP{Al) depends only 
ouAl/L. LP(Al) behaves as: (^)5/4 exp(- g^)[I - 

^ — + . . .] for asymptotically small val- 

ues of Al/L, and as {^y/^ exp{-^)[l + ci{^) + 

C2{-^y + ....] for large values. 

Moreover, we find that < A^ >= L/2. Using the 
resuh in Ref. |], Wl^ = L/12, we deduce from Eq. (2): 
< h*jy >— 5L/12, confirming the heuristic arguments 
above. Numerically the distributions of A^ and of (/i^^) 
are very similar through the whole range and have the 
same exponential decay at large values. 

We stipulate that the robustness (with respect to the 
BC) of the functional decay at large values of the MRH 
distribution is a general property of correlated surfaces. 
What is the actual functional decay for a given surface, 
and whether it is related to the distributions of other 
properties (like that of the height itself, the width, etc.), 
are left as open questions for future studies. For most 
of the surfaces we expect this decay to be exponential, 
charaterized an exponential-tail exponent a. We also sug- 
gest that a = 2 might hold for any Gaussian surface. 
Numerical simulations for the DT model |jl^ in both ID 
and 2D and for the KPZ model in 2D are clearly show- 
ing the average MRH to scale as the roughness. The tails 
of their distributions are exponential and consistent with 
the value of a = 2 . 

Having considered the late-time regime we address now 
the early-time one. In this regime the correlation length 
^ is smaller than L and growing with time. Hence the 
height variables are uncorrelated if they are separated 
by a distance Z > ^. So it is beneficial to begin by ex- 
ploring the random deposition (RD) model for which the 
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height variables are totally uncorrelated. We thus con- 
sider = H{i)-HL{t), i = l,2,...,N {N = L^) ran- 
dom uncorrelated variables. Each of them has a binomial 
distribution which, for large time, converges into a Gaus- 
sian distribution of zero mean and variance W{t) ~ i^/^. 
Extreme- values statistics then imply a scaling of the av- 
erage (or typical) MRH, to leading order, as W{t)y/2 \nN 
(the distribution of the MRH is the Gumbel distribution 
|p2[). Since the average height HL{t) is essentially un- 
correlated with each of the /i(i), the MH distribution is 
that of the MRH shifted by Ifj^. 

For realistic models with a growing ^{t), we may imag- 
ine dividing the system into TV^ = (L/^)^ cells, such 
that the surface is correlated within the cells but no cor- 
relations exist between H(r) belonging to different cells. 
Every cell has its own highest point, and they are also 
uncorrelated. Thus, we aim to apply again the extreme- 
values statistics of N^^ independent variables in order to 
estimate the MH and the MRH of the whole system. In 
this regime both have essentially the same distribution 
up to a trivial shift by the average height. This results 
from the fact that the MH of a cell has very weak (like 
a power of ^/L) correlations with the average height of 
the whole surface. 

To proceed to find their scaling behavior, we need first 
to find the functional decay of the distribution of the 
MH within every cell. To that aim, we have to deter- 
mine which of the cell average height or the cell MRH 
is more dominant {i.e. fluctuates more). For systems 
m the 2"'' cate gory their respective fluctuations may be 
estimated. The fluctuations in the average height of a 
cell are of order ^/t£,-°/^ ~ C^^', with /?' = (1 - D/z)/2. 
They have to be compared with the cell MRH fluctua- 
tions which scale as the cell roughness: ~ For 
linear models j3' = /3/2, and the MRH fluctuations are 
dominating. (3' < f3 holds for all systems of the 2"^ cate- 
gory we have checked. We conclude that for such systems 
the fluctuations in the average height may be discarded 
and therefore the decay of cell MH distribution is essen- 
tially that of the cell's MRH. 

Next we need to find the large-value decay of the 
distribution of the cell MRH. We expect it to follow 
that of P{h*^/WL) in the late-time regime, with L = ^. 
P{h*^/WL), however, is only defined for specific BC. Try- 
ing to assign such specific BC between the virtual cells is, 
of course, meaningless. Nevertheless, we will rely on the 
property, discusses above, that for all BC this distribution 
has the same Junctional decay for large values of its argu- 
ment. Therefore, this universal decay with respect to the 
BC persists even when the BC are ill-defined. It follows 
that P{h*^) of a single cell MRH is most likely to decay 
as ~ exp[{— A(/i|/W^)°}], and so does the distribution 
of the MH within every cell (recall that the cell MH and 
MRH differ by a weakly-fluctuating average height) . 

Extremal-values statistics then predict the highest of 
the maximal heights of N^^ cells to scale as: ~ 



Wj(ln Arj)!/'^. Replacing by Bt^^ and by [L/t/^/"]^ 
we obtain (after shifting by the uncorrelated HL{t)): 

<hl>=KtP[\nL-{'^)\nt + CY/'' (4) 
a 

with non- universal constants K and C. 

Two conclusions follow: (i) The growth of the RMH is 
slower than that of the roughness by a logarithmic factor. 
This is seen in the simulation results in Fig. 2 where we 
have plotted [hl/WLf vs \ogt for the ID EW model. 
The plot is consistent with the leading linear behavior 
predicted by Eq. (4). The deviation from linearity is 
in the direction expected from the next order correction 
(~ loglog ||l^) which becomes increasingly important 
near the crossover to the late-time regime (as iV^ be- 
comes smaller), (ii) The MRH depends logarithmically 
on the system size L for t < tx (when the roughness is 
L-independent). This may be observed in Fig. 3 for the 
same model. Early-time simulations for the DT model 
ID and 2D yield a similar behavior pTf . 

For systems in the first category the early-time behav- 
ior might differ. In such systems the possibility oi (3' > (3 
within a cell cannot be ruled out (in such a case the 
growth exponent of the surface assumes the value of /?'). 

This is scenario onserved in the ID KPZ model for 
which the width of the surface relative height distribu- 
tion was calculated numerically. It was found to be- 
have as exp[— (/i/W^)''], with ry = 2 at the center of the 
distribution, ?7_ ~ 1.6 in the left tail and r]+ ~ 2.4 for 
the right tail of largest heights . This behavior mim- 
ics very closely that of the velocity distribution in the 
late-time regime (with ?7_ = 1.5 and ry+ = 2.5, respec- 
tively 1^). Thus, the early-time extreme values of the 
height are determined by the fluctuations in the average 
height of a cell, and not from these of the cell MRH (in 
contrast with second category systems discussed above). 
Therefore, the power 1/a of the logarithmic term in Eq. 
(4) has to be replaced by I/77+, at least for ^ << L. 
In this regime corrections to scaling might still be im- 
portant while for larger values of ^, becomes smaller 
and the largest typical height moves towards the center 
of the distribution with possibly an effective exponent rj 
between 2.5 and 2. Similar arguments were used recently 
to explain the extremal statistics in the energetics of 
random-bond interfaces. 

In summary we have obtained the scaling behavior of 
the maximal height. In the short-time regime its growth 
differs from that of the roughness by a logarithmic cor- 
rection depending on t and L. The exponent of this cor- 
rection is determined from the exponential decay of the 
distribution of cither the velocity {V* category) , or that 
of the MRH (2"'' category), in the late-time regime. In 
this latter regime the surface is correlated throughout 
the whole system and the average MRH scales like the 
roughness. It will be most welcome if more analytical 
and numerical results for the MH and the MRH (and 
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their distributions) of other models could be obtained. 
Finally, we hope that experimental systems with accu- 
rate roughness scaling over a wide range will exhibit the 
scaling behavior of the maximal height as well. 
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FIG. 1. The distribution of HI/Wl of the ID EW model 
with periodic BC for different L. All lengths are expressed in 
terms of the lattice spacing. 



FIG. 2. {hl{t)/WL{t) f vs logt of the ID EW model for 
two values of L. The straight lines are guide for the eyes. 



Dynamics of Fractals Surfaces, edited by F. Family and 
T. Vicsek (Cambridge University Press, Cambridge, Eng- 
land, 1990). 

A. -L. Barabasi and H. E. Stanley, Fractal Concepts in 
Surfaces Growth, (Cambridge University Press, Cam- 
bridge, England, 1995). 

J. Villain and A. Pimpinelli, Physique de la Croissance 
Crystalline, (Editions Eyrolles, Paris, France, 1995). 
Dynamics of Fluctuating Interfaces and Related Phenom- 
ena, edited by D. Kim, H. Park, and B. Khang (World 
Scientific, Singapore, 1997). 

G. Foltin, K. Oerding, Z. Racz, R. L. Workman, and R. 
K. P. Zia, Phys. Rev. E 50, R639, (1994). 
Zoltan Racz and Michael Plischke, Phys. Rev. E 50, 
3530, (1994). 

M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 
56, 889 (1986). 

B. Derrida and J. L. Lebowitz, Phys. Rev. Lett. 80, 209, 
1998; B. Derrida and C. Apert, J. of Stat. Phys. 94, 1, 
1999. 

Z. Toroczkai, G. Korniss, S. Das Sarma, and R. K. P. 
Zia, Phys. Rev. E 62, 276, (2000). 

Z. Toroczkai and E. D. WiUiams, Physics Today, 24, De- 
cember (1999); J. Krug et al. P hys. Rev. E 56, 27 02, 
(1997); H. Kallabis and J. Krug, Icond-mat 980924l| ; S. 
N. Majumdar, Curr. Sci. 77, 370 (1999). 
Y. Shapir, S. Raychaudhuri, D.G. Foster, and J. Jorne, 
Phys. Rev. Lett. 84, 3029 (2000). 

E. J. Gumbel, Statistics of Extremes, (Columbia Univer- 
sity Press, NY, 1958). 

See, e.g.. Chapter I in Extreme Value Theory and Appli- 
cations, p. 15, edited by J. Galambos, J. Lechner, and E. 
Simiu (Kluwer Academic Publishers, Dordrecht, Nether- 
lands, 1994). 

S. F. Edwards and D. R. Wilkinson, Proc. Roy. Soc. Lon- 
don A 381, 17 (1982). 

F. Family and T. Vicsek, J. Phys. A18, L75 (1985). 
The same large value exponent holds for the distribu- 
tion of the range of the surface, defined as the differ- 
ence between the maximal and the minimal heights. K. 
L. Chung, Acta. Math. Hungar., 39, 65 (1982). 

S. Raychaudhuri, M. Cranston, and Y. Shapir, to be pub- 
lished. 

S. Das Sarma and P. Tamborenea, Phys. Rev. Lett. 66, 
325 (1990); D. E. Wolf and J. Villain, Europhy Lett. 13, 



FIG. 3. hi (t) vs L of the ID EW model for different values 
of t. The solid lines are fits to Ktl^ \\og L + C{t)Y'^. 
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